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Abstract 

Recent results on the entirely analytic calculation of second order corrections 
to massive quark pair production vertices induced by light quark flavours are 
reviewed. Based on the method presented in this talk the second order effects of 
the insertion of any massless one- loop vacuum polarization into the gluon line of 
the first order diagrams can be determined. The behaviour of the corrections at 
threshold for vector and axial- vector current induced massive quark production 
is discussed. 
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1 Introduction 



Within recent years a lot of work has been invested to determine two-loop radiative corrections to 
physical processes relevant for present and future collider experiments. The aim of this difficult work 
is twofold. On the one hand, one is looking for potentially large corrections and on the other one 
likes to diminish the size of the uncertainty of the corresponding one-loop calculations. Especially in 
the framework of electroweak processes involving hadrons these aspects are of importance, because 
the scale of the QCD coupling in the one-loop corrections can only be fixed by an explicit two-loop 
calculation. 

Of special interest is the process of hadron production in e + e~-coilisions. Whereas the 0(a s ) 
corrections to this process have been calculated quite a long time ago |[], |2| for all mass and energy 
assignments, the complete 0(a 2 s ) corrections have only been determined as a high energy expansion 
up to terms of 0(M A /s 2 ) Q, where yfs denotes the cm. energy and M the mass of the produced 
quarks. As far as 66-production at the Z peak is concerned this is definitely a good approximation. 
However, in view of forthcoming e + e~-collision experiments (r-charm-, i?-factory, NLC), where quark 
pairs will be produced close to the threshold, the knowledge of the entire 0{a 2 s ) corrections for all 
values of M 2 /s is mandatory. 

In our group two strategies are pursued to reach this aim. The first one consists of finding a 
numerical approximation by use of the Pade method. For that one takes information from the high 
energy regime, the threshold region and at momentum transfer zero and constructs an interpolating 
function based on conformal mapping of the physical momentum range onto the unit circle in the 
complex plane. The second strategy consists of the analytic calculation of the 0(a 2 ) corrections. 
Of course this seems to be an impossible task at the present stage. Especially certain classes of 
gluonic corrections (those which do not exist in the corresponding QED calculation) represent a great 
challenge. On the other hand, the second order corrections, which arise from the insertion of the 
vacuum polarization due to massless quarks into the gluon line of the first order diagrams (belonging 
to the so-called "non-singlet" corrections, see Fig. [l]) are accessible. In particular, the corrections 
due to massless quarks are sufficient to carry out the BLM procedure Q, which sets the scale of the 
strong coupling in the first order corrections and, if one believes in the "naive non-abelianization" 
hypothesis ||, even might provide good estimates for the complete second order corrections. 

In this talk I will report on recent developments of the analytical approach. I will concentrate 
on the calculation and the structure of the light fermionic non-singlet second order corrections. Be- 
cause the calculation is most easily performed by using dispersion integration techniques, the results 
are derived in on-shell renormalized QED. The transition to the MS scheme and to QCD is indicated 
afterwards. The presentation is kept general and can be applied to fermion anti-fermion pair pro- 
duction vertices induced by any current. In particular, I will discuss the threshold behaviour and 
the associated scales of quark pair production via the vector and the axial- vector current. I will not 
talk about massive quark pair production mediated by singlet (or "double triangle") diagrams. For 
a presentation of the Pade method the interested reader is referred to the talk by M. Steinhauser Q 
and the references given therein. 



2 Calculation of the Light Fermionic Second Order Corrections 

The non-singlet light fermionic second order corrections to the fermion anti-fermion pair production 
rate 

R FP = r^+ (-)rW+ (1) 



7T / \ 7T 



correspond to the sum of all possible cuts of the diagrams depicted in Fig. In the frame of on-shell 
renormalized QED the basic idea to calculate those diagrams is to write the one-loop corrections to 



* The reverse situation with primary production of light fermions and secondary radiation of massive ones has been 
treated in an earlier work and shall not be discussed here. Like the singlet contribution they need not to be taken 
into account for the BLM procedure as they do not contribute to the running of the coupling. 



Figure 1: Non-singlet diagrams describing primary production of massive fermion pairs (with fermion 
mass M) with additional real or virtual radiation of a light fermion anti-fermion pair (with fermion 
mass m). 



the photon propagator in terms of a dispersion relation 
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where Rff- is the Born cross-section for // production in e + e~-collisions normalized to the point 
cross-section. It is evident that eq. (||) represents a convolution of a massive vector boson propagator 
with the function Rff- The longitudinal and explicitly gauge dependent terms indicated by the dots 
are irrelevant due to current conservation. Thus the light fermionic second order corrections to Rpp 
can be written as a one-dimensional integral 
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where r^(M 2 , s, A 2 ) represents the 0(a) corrections to the FF productions rate from the real and 
virtual radiation of a vector boson with mass A at cm. energy \fs. For the special case of interest, 
m — > 0, the evaluation of the integral (|3|) can be simplified enormously by subtracting and adding Rff 



at its asymptotic high energy value R^ 
rf(M 2 ,s) = I 
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Because Rff reaches its asymptotic value already for A 2 much smaller than the hard scales M 2 and s 
we can replace r^-\M 2 , s, A 2 ) in the first integrand on the r.h.s. of eq. @ by A l \M 2 , s, 0), the one- 
loop contribution to Rpp in eq. (|T|). The Bloch-Nordsieck theorem Q assures that the corresponding 
limit exists. The first integral can then be trivially expressed in terms of the moment 
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The evaluation of the second integral constitutes the main effort. The final result for the second order 
corrections due to one light fermion / can be written in the form (x = f) 
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where 5^ is a complicated function of the ratio M 2 /s involving Tri- (Lis), Di- (L12) and usual 
logarithms 



1C]. For the light fermionic corrections we are interested in the moments read 

RL = i. r/ 
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The result in the on-shell scheme can now be easily transferred to the MS scheme by expressing the 
fine structure constant a in terms of the MS coupling 
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which effectively replaces ln(m 2 /s) in (0) by ln(/u 2 /s). At this point I would like to emphasize that 
the moments can also be extracted directly from the vacuum polarization function due to a massless 
fermion anti-fermion pair (x = /), 
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To arrive at the corresponding light fermionic second order corrections in the frame of QCD we 
have to multiply the corresponding SU(3) group theoretical factors N c = 3, T = 1/2 and Cp = 4/3. 
In particular the moments for light quarks in QCD read -R^ ,QCD = T R^ (n = 00, 0). 



3 The Method of Moments 

The method described above can be applied to determine the corrections from the insertion of any 
massless vacuum polarization into the gluon line, as long as the absorptive part of the vacuum polariza- 
tion function approaches a constant in the high energy limit. This is equivalent to the occurrence of at 
most one single power of the logarithm ln(— q 2 / n 2 ) in the renormalized vacuum polarization function 
in D = 4 dimensions and is certainly true for any one-loop vacuum polarization function. Thus we can 
easily determine for example the second order corrections from the insertion of a vacuum polarization 
due to massless coloured scalars into the gluon line. According to eq. @ the scalar moments read 

« CD = T Q) , ^ QCD = T Q ln4 - |) . (10) 

We are also in a position to calculate the second order gluonic self-energy contribution to massive 
quark production, as illustrated in Fig. ||, by determining the corresponding moments of the gluonic 
and ghost contributions to the 0{a s ) gluon propagator corrections, 

The gauge parameter £ is defined via the gluon propagator in lowest order 

i(- r+i tf\. (12) 

Evidently the gluonic moments and r^' QCI> , defined in analogy to eq. (^), are not gauge invariant. 
However, for £ = 4 the combination a 2 R^ CD /(Sir) coincides with the complete gluonic 0(a 2 ) contri- 

( 2) OCD 

bution to the QCD /^-function, and thus r g \(,=i accounts for the leading logarithmic behaviour of 
the sum of all gluonic 0(a 2 ) diagrams in the high energy limit. It is not surprising that such a gauge 
can be found. However, it is a remarkable fact that for the same choice of £ also the gluonic contri- 
butions to the perturbative QCD potential [|ll]] can solely be expressed in terms of the moments (|ll|). 
This point will be examined in more detail in the next section. Further application of the method of 
moments to corrections to Rpp beyond the second order in the strong coupling can be found in g, 




4 The Threshold Behaviour 



The examination of the threshold region is of special interest. Here, the contributions from the 
exchange of longitudinal and transversal polarized gluons can be clearly distinguished. The former 
gives rise to the instantaneous Coulomb interaction which cancels one power of (3 = \/l — 4M 2 /s 
from the phase space and leads to the well-known finiteness of the vector current cross-section at 
threshold. The latter is responsible for the relativistic corrections. Both contributions are governed 
by different momentum scales. Prom general considerations one can infer that the scale relevant for the 
instantaneous Coulomb interaction should be of the order of the relative momentum of the produced 
massive quarks, y/s/3, whereas the scale for the hard/transverse corrections should be of the order of 
the cm. energy, y/s. However, only an exact calculation of the non-singlet light fermionic second order 
corrections allows for an accurate determination of those scales via the BLM scale setting procedure. 

Taking in account only the second order effects of rif light quark species the threshold expansion 
of the massive quark pair production rate induced by the vector current reads 



R v = N c C F -a s (^)7r 



1+ sl * 



1 fa s (n 
3 



N C C F -asiji 2 )^ 



RL QCD In 



16 fi 



2 + 4 )+^o 



/,QCD 



1 U f^massless( a e ^^) 



P 



P 



(13) 



including terms up to 0(/3). In the ^-scheme, where the strong coupling is defined as the all order 
effective charge in the QCD potential 
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the corresponding BLM scales can be directly read off the arguments of the vacuum polarization 
functions Tl^^ sless on the r.h.s. of eq. (13). The reader should note that the second line of eq. (|i~5|) 
is a new observation || which has not been made in the original paper [11]. Obviously the leading 
second order threshold contribution to R v from massless quarks is in one-to-one correspondence to 



the form of the QCD potential. This is in accordance with non-relativistic quantum mechanics. As a 
consequence the leading second order gluonic threshold contributions to R v , which are proportional 
to the colour factor Ca (and represent purely non-abelian contributions) can be inferred direcly by 
adding n^ 1 ^ z (s /3 2 )|^=4 on the r.h.s. of eq. (13). A similar conclusion is impossible for the hard 



corrections, because they represent effects which cannot be described by the instantaneous Coulomb 
potential. 

Taking into account only the second order non-singlet effects of rif light quarks flavours the 
axial-vector induced massive quark pair production rate at threshold can be cast into form 
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The reader should note that the leading contribution to R is of 0((3 2 ). In analogy to the vector 
current case the BLM scales in the V^-scheme can be identified from eq. (|16|). The 0((5 2 ) second order 
gluonic corrections oc Ca can be determined in complete analogy to the vector current case. A formal 
proof of this statement will be presented in later work. 
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